Effects of chiral restoration on the behaviour of the Polyakov loop at strong coupHng 
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We discuss the relation between the Polyakov loop and the chiral order parameter at finite tem- 
perature. For that purpose we analyse an effective model proposed by Gocksch and Ogilvie, which 
is constructed by the double expansion of strong coupling and large dimensionality. We make im- 
provements in dealing with the model and then obtain plausible results for the behaviours of both 
the Polyakov loop and the chiral scalar condensate. The pseudo-critical temperature read from the 
Polyakov loop turns out to coincide exactly with that read from the chiral scalar condensate. Within 
the model study based on the strong coupling expansion, the coincidence of the pseudo-critical tem- 
peratures results from the fact that the jump of the Polyakov loop in the presence of dynamical 
. . . quarks should signify the chiral phase transition rather than the deconfinement transition. 
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Introduction Quantum Chromodynamics (QCD) is commonly accepted as the fundamental theory of the strong 
Q ] interaction. It has been intensely argued how a thermodynamic system described by QCD goes through the phase 
transitions at sufficiently high temperature, namely, the colour deconfinement transition and the chiral phase transition 
around the critical temperature ^150 McV [|l[. 
, The analytic study of the deconfinement transition was first offered within the framework of the strong coupling 
■ expansion on the lattice |^ . In the absence of dynamical quarks a pure gluonic system has the centre symmetry and 
', the Polyakov loop provides a criterion of confinement . For the purpose of looking into the spontaneous breaking of 
^ ' the centre symmetry, the effective action in terms of the Polyakov loop variables has been constructed by the strong 
\ coupling expansion, as well as by the perturbative calculation. It is found that the resultant effective action leads to a 
I • second-order phase transition |^ for the SU(2) gauge theory and a first-order phase transition for the SU(A^ > 3) 
CO ' gauge theories. These are consistent with the anticipation from the point of view of the universality classifications 
, as well as of the lattice observations |^ . Also from the analysis of the effective action it has been shown that the 
perturbative vacuum of the pure gluonic system becomes unstable beyond some critical coupling strength . 

When dynamical quarks in the fundamental representation are present in a theory, the centre symmetry is broken 
explicitly and the Polyakov loop is no longer regarded as an order parameter to characterise the deconfinement 
transition Then, the chiral symmetry associated with light quarks plays an important role in hadronic properties. 
There are many effective approaches based on the chiral symmetry, such as the linear-sigma model, the Nambu-Jona- 
O i' Lasinio model, the chiral random matrix model and so on In those model studies, features of colour confinement 
D , are hardly taken into account. The relation between the confinement and the chiral symmetry breaking has been 
' discussed theoretically in a variety of contexts such as the hclicity conservation |^ , the anomaly matching condition pX| ] 
^ , and so on JTTI . The difficulty in studying their interrelation lies in the fact that no proper criterion of confinement is 
t'"! ' established so far, despite of great efforts |p^ . 

rS [ In the meanwhile, it has been found in the lattice simulations that the Polyakov loop behaves like an approximate 
order parameter even in the presence of dynamical quarks. Then the jumps of the Polyakov loop and the chiral 

" " ' condensate are observed at the same critical point [|l3j . Precise analyses thereafter in which the pseudo-critical points 
are defined by respective susceptibilities have sustained this coincidence |Q . It is often said that the deconfinement 
and the chiral restoration should take place simultaneously, though the physical reason for the coincidence is still 
obscure. 

The author of Rcf. [ p5| proposed the following simple explanation for the simultaneous jumps; at low temperature 
where the chiral symmetry is spontaneously broken, the explicit breaking of the centre symmetry is suppressed by 
the heavier mass of constituent quarks rather than the lighter mass of current quarks. Consequently the expectation 
value of the Polyakov loop stays small in the confined or chiral broken phase at low temperature. Once the constituent 
quark mass drops off in the chiral symmetric phase at high temperature, the expectation value of the Polyakov loop 
no longer receives such suppression. It follows that the jump of the Polyakov loop signifies not the deconfinement 
transition but the chiral phase transition. The simultaneous jumps are observed simply because the behaviours of 
both the Polyakov loop and the chiral condensate indicate a single phase transition of the chiral restoration. 

This scenario had been partially confirmed within an effective model constructed by the double expansion of strong 
coupling and large dimensionality on the lattice m&. In Ref. ffffl Gocksch and Ogilvie found that the given model 
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leads to the deconfinement transition at Td = 175 MeV for the pure gluonic system and the chiral restoration at 
T-^ = 471 MeV. The authors also argued that the constituent quark mass tends to suppress the magnitude of the 
Polyakov loop, though the simultaneous jumps could not be reproduced there. Actually in Gocksch and Ogilvie's 
results, the transition temperatures relevant to the deconfinement and the chiral restoration are too different to affect 
each others. 

The purpose of this letter is to show that the jump of the Polyakov loop should be attributed to the chiral phase 
transition rather than the deconfinement transition. To that end, we reexamine the model given by Ref. ]l6| , which 
we call the Gocksch-Ogilvie model hereafter. We make improvements in phenomenological approximations employed 
in the original work and then find the deconfinement transition at Td = 208 MeV for the pure gluonic system and 
the chiral phase transition at = 162 MeV in the absence of the Polyakov loop dynamics. As compared with 
the empirical values these transition temperatures are more reasonable than those originally derived in Ref. [ p^ . 
Furthermore, we analyse the coupled dynamics of the chiral order parameter and the Polyakov loop in a numerical way 
and demonstrate that the pseudo-critical temperature in regard to the Polyakov loop coincides with the pseudo-critical 
temperature of the chiral restoration at Tc — 187 MeV. 

Gocksch-Ogilvie model The effective action of the Gocksch-Ogilvie model is constructed by the double expan- 
sion of strong coupling, where the confining property is almost trivial, and large dimensionality, where the dynamical 
chiral symmetry breaking is realised on the lattice ||l^ . Here we present only the expressions from Ref. ; 

Sos[L,X] = -J^TrcL(n)Trcit(m) + ^ ^ A(n)F(n, m)A(m) 



\T.{^MooMNrE)] +Tr.ln[l-H ^^^^^|^^^^ (X + Lt)]}, (1) 



where X{n) is the meson field, the condensation of which is responsible for the dynamical breaking of the chiral 
symmetry. Since the model is based on the staggered formalism on the lattice, the flavour contents of meson fields 
are resolved in the differing lattice sites. A^c is the number of colours. L{n) denotes the Polyakov loop defined in the 
d-dimensional space-time by 

L{n) - Yl Udin,n4) (2) 

n4 — 1 

with Nr lattice sites in the thermal (temporal) direction. The Polyakov loop is a matrix in the colour-space and Trc 
represents the trace over the colour indices. The static energy E of quasi-quarks is given by 



E = smh-UJ- — -X + m] (3) 



with the current quark mass denoted by m. In the above expression {d— 1) is the number of spatial dimensions. The 
propagator V{n, m) of meson fields is given by 

""^ " 2(d- 1) ^ + ^n,m~^ ' (4) 

where j runs over only the spatial directions. The strength of the nearest neighbour interaction, J, is determined as 
a function of the temperature, in other words, the temperature is specified by J: In the strong coupling limit (J <JC 1) 
we can express the correlation function of the Polyakov loops in twofold ways as 

(Trcit(ri)Trci(m)) - jl™""! ~ Q-li'yo.\in-n\ ^ 

where a is the string tension and a is the lattice spacing. Accordingly the temperature in the Gocksch-Ogilvie model 
is fixed as 

1 era 

It should be emphasised that the coefficient 1/2 in front of the last part of Eq. (|l|) was missing in the original 
work by Gocksch and Ogilvie. This factor is important in counting the number of flavours from the phenomenological 
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point of view. But for the coefficient 1/2, the scalar field A(n) generates too many mesons with four degenerate light 
flavours peculiar to the staggered formalism. This coefficient effectively reduces the model to that with only two {u 
and d) light quarks, as often adopted in lattice simulations. 

The prescription of taking the square root of the Dirac determinant is reliable only in the weak coupling limit where 
the Dirac operator is well localised. In the present case, the prescription is also reliable, though the model is based 
on the strong coupling expansion. The counterpart of the Dirac operator is absolutely local because the last part of 
Eq. (|l|) comes from the integration with respect to the quark field without its kinetic term which is absorbed in that 
of mesons. 

With the factor 1/2, the formulae at zero temperature given in Ref. [|l^ should receive modifications. Certainly the 
expectation value of A at zero temperature becomes •y/l/2 from 1, but all the physical quantities such as the hadron 

masses and the chiral condensate are not affected if just the current quark mass is divided by Actually there 

are two parameters inherent in the model, namely, the lattice spacing a and the current quark mass m. According to 
the formulae given in Ref. with modifications from the factor 1/2, the model parameters can be fixed as follows 
so as to reproduce the pion mass 140 MeV and the p meson mass 770 MeV: 

a^^ = 432 MeV, m = 5.7 MeV. (7) 

We would emphasise the following point: In the zero temperature analysis like Ref. jlTt , too many light flavours 
would cause no serious problem as long as the mass spectrum in only the u and d quark sector is concerned. In the 
chiral limit (m — 0) the mean-field approximation leads to the same answer for each quark flavour regardless of the 
flavour number. In the presence of small current quark mass, the difference of the flavour number can be absorbed in 
the difference of the current quark mass. At flnitc temperature, on the other hand, how many degrees of freedom are 
excited at that temperature should affect signiflcantly the relevant temperature to the chiral restoration. In this case, 
since the number of pions (Nambu-Goldstone bosons) is too large, the resultant potential energy prevents the chiral 
restoration until considerably high temperature is reached.^ Thus the factor 1/2 is important in order to estimate 
the chiral transition temperature quantitatively. 

Mean- field approximation Following the treatment of Ref. [|l6j , we adopt the mean-field approximation to deal 
with the effective action (|l|). As to the meson field, we simply replace it by a condensate A. The Polyakov loop L{n) 
is integrated out with the mean-field action [ p^ , 

Smt[L] - -| ^{Trei(n) + Tv.L^n)}, (8) 

n 

where x is the variational parameter to be determined afterwards from the extremal condition on the free energy. The 
mean-field free energy is then given by 

/3/mf(x,A) = -Ar-('^-i'((-S'eff[L,A] -f 5„,f[L])„,f + ln y 2?Le-^»^[^l) 

= /3/if(^)+/3C(2;,A), (9) 

where (• • • )mf stands for the expectation value calculated by means of the mean-field action (||). In the expression of 
the mean-field free energy, the physical implication of each term is plainly understood as follows: In the first line, the 
first term is the internal energy tending to make the system ordered and the remaining terms can be regarded as the 
entropy tending to make the system disordered, f^fix) represents the pure gluonic part of the mean-field free energy 
and f^fix, A) corresponds to the chiral part with the coupling between the Polyakov loop and the scalar condensate. 
They are explicitly written as 

(3fUx) = -2(d " In Hx)) ' + In H^)) , (10) 

^.a / 7x NcNr^n iVc , r w M 1 i(x:cosh(NrE)) , , 

/3/mf(^,A) = ^^A^ - ^ln[cosh(iV.£;)] - - ^ (11) 



^ This might seem contradict to the fact that the chiral transition temperature is lowered as the number of flavours increases. In the 
treatment of the Gocksch-Ogilvie model, however, we exploit the mean-field approximation and any fluctuation of pions, which is 
responsible for reducing the transition temperature, is frozen from the beginning. Thus the potential energy alone is involved in the 
analysis and it gives higher transition temperature with more flavours. 
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where I{x) is defined by 

lix) = fdL exp[|Trc(L + it)] ^ Arj ^ det (x), (12) 

m 

and I(x; a) is given by 

i{x;a) — ^iii2---ijvc^iii2-- jNc^™-H+ji (^) • ■ ■ -^™-«a+ia (2^! Q^) • ■ ■ 

a— 1 m— — 00 

/„(x;a)= r"^ln[l + ^je-™^*+-^. (13) 
./n 27r L a J 



Here /„(a;) denotes the modified Bessel function of the first kind. Once we expand /„(x;a) in terms of 1/a and 
furthermore expand the above expressions in terms of x, we can immediately reproduce the results of Ref. JlGf . These 
expansions in terms of 1/a and x are not always regarded as reliable in fact. Although the expansion in terms of 
X works well for the estimate of I{x), it is not a good approximation for f^i{x) due to significant cancellations in 
the right-hand-side of Eq. ([lo|). The expansion in terms of 1/a becomes uncertain when the expectation value of the 
Polyakov loop is large. Therefore we do not adopt such expansions in examining the mean-field free energies in the 
present work. 

The deconfinement transition in the pure gluonic system is described by f^{{x), while /j„f (x, A) determines the fate 
of the chiral symmetry in the presence of the Polyakov loop. Actually f^[{x) results in the first-order phase transition 
at 2{d— l)Jc = 0.806 in the case of A^c = 3. The corresponding temperature read from Eq. (^ is Td = 208 MeV with 
the numerical value of the string tension a = (425 MeV)^ substituted. 

Then, we will look into the chiral dynamics somewhat closely in an analytic way in the case of the chiral limit 
(to = 0). When the phase transition is second-order, that is the case in the present analysis, the transition temperature 
can be deduced from the condition that the quadratic term with respect to the order parameter changes its sign, or 
passes across zero. Thus we can derive the analytic expression for from the expansion of f^i{x, A) in terms of A. 

The calculation is simplified when the Polyakov loop is put as trivial, i.e., i ~ 1 for the moment. Then we can 
approximate the last term of Eq. ( pi] ) by —{Nc/2) ln[l + 1/ cosh(A^T-£')] looking back at Eq. (^, which leads to 

From this expansion, the chiral phase transition temperature in the absence of the Polyakov loop dynamics is imme- 
diately determined as 

= ^^—^a-^ = 162 MeV (15) 
8 

with the numerical value of a given by Eq. (0). The above expression is essentially the same as (3.23) in Ref. ||l^. 

Another interesting limit in which the calculation is simplified is that the Polyakov loop is forced to vanish, i.e., 
X = 0. Then it follows that 

(Tr. ln[l j^J^^iL + ^t)] )^^^^ . -A^, ln2 + Vl^A^.A^.|A| - l^Llil^p. (16) 

The quadratic term proportional to A^ is exactly cancelled by the contribution from Ac ln[cosh(A^7.i?)] ~ {d — 
l)AcA^A^ /4 and only the tree-term remains. What is surprising is that the linear term of A appears 

in Eq. dlq). Owing to the presence of this linear term, the stationary point with respect to A always leaves from 
zero. In other words, the chiral symmetry is broken at any temperature. If the vanishing Polyakov loop, x ~ 0, 
has something to do with confinement even in the presence of dynamical quarks, this result suggests that the chiral 
symmetry must be broken in the confined vacuum, which is in agreement with the arguments in Refs. |l0| . 

Hence, the chiral dynamics is substantially affected by the Polyakov loop dynamics. It means at the same time 
that the behaviour of the Polyakov loop is deeply related to that of the chiral order parameter. 

Numerical results We must search for the global minimum of the free energy fmf{x. A) numerically in general, 
except for the above special cases of to = 0, L ~ 1 and to = 0, x = where analytic treatments are feasible. The 
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variational parameter, a;, and the scalar condensate, A, are determined as functions of temperature. Then we can 
readily acquire the expectation value of the Polyakov loop by using (see Eq. (|l^)) 




The results are shown in Fig. ^ It is apparent that the behaviours like phase transitions are observed around the 
temperature ~ 180 MeV in Fig. ||. As compared with the pure gluonic result = 208 MeV (shown by the dotted 
curve for reference) it seems that the critical point concerning the Polyakov loop dynamics is smeared by dynamical 
quarks and, as a result, two crossovers in terms of the Polyakov loop and the scalar condensate are observed nearly 
at the same point. 

To make the argument more quantitative, we can define the pseudo-critical temperatures by means of the peak of 
susceptibilities in a similar way to Ref. p^ . Here we make use of the simplest ones, which are immediately feasible in 
the present analysis, that is, the temperature susceptibilities Xt = d{TrcL)/dT/Nc and Xt — —dX/dT. The results 
are shown in Fig. || The pseudo-critical temperatures are found exactly at the same point Tc = 187 MeV. As to 
the moderate peak of Xt^ slightly above ~ 200 MeV, it should be regarded as an accidental artifact because this peak 
would vanish if the current quark mass m is raised or lowered by hand. 

The coincidence of the pseudo-critical temperatures completely agrees with the lattice observations. Within the 
Gocksch-Ogilvie model, the jump of the Polyakov loop does not mean the deconfinement transition. Actually it is 
caused by the last term of Eq. (|l|) through which the behaviour of the Polyakov loop would reflect that of the scalar 
condensate. In analogy with the hopping parameter expansion (see Eq. (2.11) in Ref. we can regard 1/ cosh(A'^r£') 
as the strength of an external field to break the centre symmetry. The behaviour of the Polyakov loop is governed 
effectively by the chiral dynamics through this strength of an external field. 

We consider that the exact coincidence of the pseudo-critical temperatures in the present study provides a credible 
support for the argument of Ref. : The behaviour of the Polyakov loop indicates only the chiral phase transition 
rather than the deconfinement transition. Thus the coincidence of the critical temperatures is only a consequence of 
viewing a single phenomenon, i.e., the chiral restoration. 

Finally we shall comment upon the plans for further work. Since the number of the space-time dimensions is four at 
most, it is necessary to make sure the convergence of the large dimensional expansion and to confirm the quantitative 
success in describing the Polyakov loop and the chiral dynamics even with next-to-leading order contributions taken 
into account. 

It would be intriguing to construct a similar effective model in the formulation of the Wilson fcrmion [2l|| . Then 
the flavour decomposition becomes straightforward at the cost of the intricate structure of the Dirac indices. Also an 
extension to introduce adjoint fermions, which does not break the centre symmetry, would be interesting. In contrast 
to the case of fundamental fermions, two distinct transitions are found in the lattice simulations p^ . It would be 
challenging to reproduce such results within the model study parallel to the present treatment. Such analyses would 
serve as a double-check for the speculation on the behaviour of the Polyakov loop discussed in this letter. 
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FIG. 1: The behaviours of the order parameters as functions of the temperature. Transition-like jumps are observed simulta- 
neously around Tc ~ 180 MeV. The dotted curve is the result in the pure gluonic case for reference. 



FIG. 2: The temperature susceptibilities for the Polyakov loop and the chiral scalar condensate. 
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